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Abstract 
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I. INTRODUCTION 



The next decade promises to be a very exciting time for strong interaction physics. 
With ever increasing computer power and impressive progress in developing new techniques 
to simulate quantum field theories, one hopes that lattice simulations of simple hadronic 
systems will provide rigorous and reliable predictions of QCD for strong interaction observ- 
ables. While one looks forward to fully- unquenched QCD simulations performed with the 
physical values of the light-quark masses, rrig, such simulations are presently prohibitively 
time-consuming. At present, and for the foreseeable future, unphysical theories will be simu- 
lated because, in contrast with QCD, the simulations can be performed in a reasonable time 
frame A second motivation for simulating unphysical theories is particular to nuclear 
physics, and is related to the unnaturally large values of the S-wave scattering lengths. We 
will discuss this point in the final section of the paper. A commonly simulated unphysi- 
cal theory is quenched QCD (QQCD) where disconnected quark-loop diagrams (the quark 
determinant) are omitted. While QQCD simulations of strong interaction observables can 
be performed with small rriq, they have the distinct disadvantage of not being related to 
QCD except in the large- A^^c limit 0. A more interesting unphysical theory is partially- 
quenched QCD (PQQCD) in which the quark masses, ms (the "S"' stands for 
"sea"), used in evaluating the disconnected quark- loops are larger than the masses of the 
quarks connected to external sources, my (the "V" stands for "valence"). By computing 
strong interaction observables in PQQCD and performing the extrapolation ms,mv rriq 
one recovers the QCD observables that one is interested in. It is for this last step — the 
extrapolation to the physical values of rriq — that effective field theory (EFT) is required. 
The EFT's describing QCD in the low-momentum regime in the pseudo-Goldstone boson 
sector (chiral perturbation theory, xPT), and in the single baryon sector (heavy baryon 
chiral perturbation theory, HBxPT), are well established, and their extension to PQQCD 
in the form of PQxPT P,^,|^,PJ^ and PQHBxPT []8|,|9|JlO[| have been accomplished relatively 
recently. 

The construction of an EFT to describe the low-momentum dynamics of multi-nucleon 
systems has proven to be extremely challenging. In the very low-momentum regime, where 
the typical momentum of the external particles involved in a given process is much less than 
the mass of the pion, p <^ m^r, and hadronic production is therefore kinematically forbid- 
den, an EFT, EFT(7f) , can be constructed from nucleons and photons (and any 
other low-momentum transfer probes) quite simply. The fact that there is a bound state near 
threshold in the ^Si — ^Di coupled- channels, and a pole on the second-sheet near threshold 
in the ^5*0 channel means that at least one operator in the EFT(7f) Lagrange density must be 
treated non-perturbatively in these channels. The choice of operators to be resummed and 
details of the perturbative expansion are defined by the power-counting in EFT(7f). Despite 
chiral symmetry not being a good symmetry for p -C rriT^, isospin remains a good symmetry. 
The only input into the construction of EFT(7f) is Lorentz invariance, electromagnetic gauge 
invariance, baryon number conservation and the approximate isospin symmetry, the break- 
ing of which can be included perturbatively. In the kinematic regime where the momenta 
involved in a given process are larger than m^, the pion must be included as a dynamical 
field. It was Weinberg's pioneering efforts |]I5| in the early 1990's in this kinematic regime 
that initiated interest in developing EFT for nuclear physics. Weinberg attempted to con- 
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struct an EFT for nuclear processes and nuclei involving momenta all the way up to the 
chiral symmetry breaking scale A^, and necessarily included the pion as a dynamical degree 
of freedom. The power-counting that he developed, known as Weinberg power-counting 
(W), involves a chiral expansion of the nucleon-nucleon potential using the same power- 
counting rules that are used in the meson and single nucleon sectors. The chirally expanded 
potential is inserted into the Schrodinger equation to determine observables, such as phase 
shifts. Unfortunately, there is a formal problem with this power counting in some of 
the scattering channels, particularly the ^5*0 channel. However, extensive phenomenological 
studies with W power-counting appear to be in good agreement with data |r7| , p!^ , pD| , pD| , 
where such comparisons are possible, and the formal problems appear to have little impact 
when a massive regulator is used with a mass scale that is not radically different from a 
few hundred MeV. The formal problems with W power-counting led Kaplan, Savage and 
Wise (KSW) to develop a power-counting in which the momentum-independent four- 
nucleon operator is promoted to one lower order in the chiral expansion, and consequently 
pion exchanges are subleading and treated in perturbation theory. This power-counting is 
formally consistent and gives renormalization group invariant amplitudes order-by-order in 
the EFT expansion. However, Fleming, Mehen and Stewart (FMS) showed that the 
scattering amplitude in the ^5*1 — ^Di coupled channels diverges at next-to-next-to-leading 
order (NNLO) at relatively small momenta and KSW power-counting fails. FMS found that 
a contribution that remains large in the chiral limit destroys the convergence: it is the chiral 
limit of the tensor force that "does the damage" . Recently, it was suggested that one should 
expand observables about the chiral limit p3| (BBSvK power-counting). BBSvK power- 
counting has all the nice features of W and KSW counting: the chiral limit of the tensor 
force is resummed at leading order (LO) along with the momentum- and mg-independent 
four-nucleon operator in the ^5*1 — ^Di coupled channels, while pions are perturbative in the 
^5*0 channel, and in higher partial waves, where analytic calculations are possible. 

In recent work we showed that hairpin diagrams in PQQCD give rise to a component 
of the nucleon-nucleon (NN) potential that falls exponentially at long-distance and therefore 
does not have the Yukawa behavior found in QCD. Thus, measuring the long-distance be- 
havior of the NN potential p5|j26|j27|j28[] does not provide information about QCD unless the 
hairpin contribution can be removed in a rigorous way. While the presence of this behavior 
is quite discouraging, one might focus on the behavior of S-matrix elements rather than on 
the NN potential itself (for a recent survey of the status of lattice calculations of the NN 
potential see Ref. [^). In this work we develop the partially-quenched EFT that describes 
the two-nucleon sector using BBSvK power-counting. In channels with higher partial waves, 
we give analytic expressions for scattering amplitudes and a few characteristic scattering 
volumes to next-to-leading order (NLO) in the partially-quenched EFT. In the ^^o-channel 
we give analytic expressions for the scattering amplitude, scattering length and effective 
range to NLO in terms of the valence and sea quark masses. As in QCD, the ^5*1 — ^Di 
coupled channels in the partially-quenched EFT are somewhat more complicated; we pro- 
vide the NN potential at NLO that is required to generate the NN phase-shifts, 5o,2 and 
mixing-parameter, ei, by solving the Schrodinger equation. 
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II. THE PARTIALLY-QUENCHED EFT CALCULATION 



In BBSvK power-counting S-matrix elements are an expansion about the chiral limit, 
where the expansion parameter, Q, is Q ~ 1/3 ~ tjIt^/Ann, where the constant A^vat = 
87if'^/{g\Mfyf) is determined by the relative size of pion exchange, and m.,^ is the physical 
value of the pion mass. In the two S-wave channels, the momentum and mg-independent 
four-nucleon operators, along with the chiral limit of one-pseudo-Goldstone-Boson-exchange 
(OPGBE), are resummed to all orders (each iteration is the same order in Q) and this 
sum constitutes the LO scattering amplitude. In the ^S'o-channel this corresponds to KSW 
power-counting and the scattering amplitude can be computed analytically, while in the 
^5*1 — ^Di coupled channels, the scattering amplitude must be determined numerically. In 
the higher partial waves, the chiral limit of OPGBE provides the LO contribution to the 
scattering amplitude as the four-nucleon operators are suppressed by additional powers of 
momentum. 

It is the ultra-violet (UV) behavior of the theory that requires an expansion about 
the chiral limit, and in particular allows control of the very singular diagrams as r — > 
in coordinate-space. The deviations from the chiral limit in OPGBE that formally occur 
at higher orders in BBSvK power-counting are UV safe and can therefore be included at 
NLO without compromising the renormalizability of the theory. The price for not including 
them is that the long-distance behavior of the theory must be recovered order-by-order in 
perturbation theory and convergence is somewhat slow [^. By contrast, the two-pseudo- 
Goldstone-Boson-exchange (TPGBE) diagrams are singular away from the chiral limit, and 
therefore only the chiral limit can be retained at NLO; keeping the full TPGBE introduces 
divergences that cannot be renormalized at NLO. In this work we compute to NLO in the 
EFT with BBSvK power-counting. In the S-wave channels at NLO there are contributions 
from OPGBE (the full meson mass dependence is retained), from momentum- and rriq- 
independent four-nucleon operators, from the leading momentum-dependent four-nucleon 
operators (p^) and from the four-nucleon operators with a single insertion of mg. In the 
higher partial waves, the four-nucleon operators contribute beyond NLO and thus only 
OPGBE contributes at NLO. 

We will work in the isospin limit of the SU{4:\2)l ® SU{4\2)ii PQ^PT. This means that 
the u, d, u and d quarks in the valence and ghost sectors are degenerate and the j and 
I quarks in the sea sector are degenerate. The formalism for this theory can be found in 
Refs. and we will not describe it here. 



A. Partial Waves with L > 



At NLO (0(Q°)), the partial waves with L > (higher partial waves) receive contribu- 
tions only from OPGBE, as shown in Fig. [I|, and no resummation of diagrams is required. 



In PQxPT the potential between two-nucleons due to OPGBE is ||2^ 
arising from the interaction Lagrange density in PQQCD (in the isospin limit) 
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The TT and rj propagators are of the form 
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where / ~ 132 MeV, is the isovector axial couphng constant and 5^0 is the isoscalar axial 
coupling. The mass mss is that of a meson composed of two sea quarks while m^r is the 
pion mass which is, of course, composed of two valence quarks. The propagators in eq. @ 
clearly exhibit the correct behavior in the QCD limit, where the coefficient of the double- 
pole contribution in the rj propagator vanishes, and the single pole contribution is absent. 
By treating m^r, mss and |q| all of 0{Q), both OPGBE contributions are the same order in 
the power-counting. 



(a) 



(b) 



FIG. 1. The LO contribution, 0{Q^), to scattering in the higher partial waves from OPGBE. 
Diagram (a) corresponds to the exchange of vr, while diagram (b) corresponds to the exchange of 
■q with a double-pole propagator, as given in eq. (^), denoted by "X ". 

In QCD, the scattering amplitudes in the higher partial waves, Aj^^^^si total angular 
momentum J, total spin 5*, and initial and final state orbital angular momentum L and L' 
respectively, are well known to 0{Q) [0. In the spin-singlet channel {S = 0) the scattering 
amplitude in partial waves with J = L > 0, Aj^ff^j, is 
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Ref. g. 



[z) is an irregular Legendre Polynomial of order n, using the conventions of 
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and where the variable z is given in terms of the pion mass and nucleon center-of-mass 
momentum as 2; = 1 -|- rn^/{2p'^). The quantity / is the isospin of the channel under consid- 
eration. In the spin-triplet channel {S = 1) the expression is somewhat more complicated 
due to the fact that the spin and orbital angular momentum can couple to produce three 
different total angular momentum states, J = L — 1, L, L + 1. It is straightforward to show 
that the amplitudes for scattering between states of the same orbital angular momentum 
are 
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where the constant 8(2^ is given by 
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The amphtudes for scattering between states with orbital angular momenta that differ by 
two units, AL = 2, induced by the tensor component of the interaction are 

4SS2,i = ^ ^/2"^' [ QlM^) + QLiz) - 2Ql+,{z) ] , (8) 

where 8(2^^'^ is given by 

- 2J+1 ■ 

In order to arrive at the partially-quenched amplitudes it is convenient to note that the 
contribution from single r) exchange can be obtained from OPGBE by taking a derivative of 
with respect to m^, and multiplying by appropriate constant factors. Generically, 

^ffiB.= (lMH'(2/ + l)^(^^%^)|)^S?^']. . (10) 

It is straightforward to show that the partially-quenched amplitudes in the spin-singlet 
higher partial waves are 
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and in the spin-triplet higher partial waves are 
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The phase-shifts in the spin-singlet channels can be easily extracted from the scattering 
amplitudes given in eq. ([TTl) by using the relation (for non-relativistic systems) 



^j,j.jfl = TT log + i-7r^^JJJ,o ) ) (13) 



2i " V 27r 

from which parameters in the effective range expansion can be determined for p < m^, mss- 
It is somewhat more complicated to determine the phase-shifts in the spin-triplet channels 
as one has to disentangle them from the mixing parameters. However, at NLO {0{Q^)), 
the mixing effects are higher order in the EFT and one can straightforwardly determine 
parameters in the effective range expansion. In the P-waves, the scattering volumes, defined 
to be 

a( Pj) = - \im , (14) 



are found to be, at NLO, 
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for which the QCD limit agrees with the well-known results 30 



B. The % Channel 

The scattering amplitude in the ^Sq channel can be determined analytically order-by- 
order in perturbation theory as BBSvK power-counting coincides with KSW power-counting 
in this channel. The momentum and mg-independent four nucleon operator with coefficient 

Cq "^"^ enters at LO, and the bubble chains that it generates, as shown in Fig. ^ are resummed 
to all orders to produce the LO scattering amplitude [^]. At NLO there are several different 




FIG. 2. The LO contribution, 0{Q ^), to the scattering amplitude in the ^Sq channel. 

contributions. There is a contribution from OPGBE that can be dressed in a variety of 
ways by the LO amplitude as shown in Fig. |^, and each dressing remains 0{Q^). There 
is a contribution from a momentum-dependent (p^) operator with coefficient C2 that is 
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FIG. 3. The NLO contributions, 0{Q^), to the scattering amplitude in the ^5*0 channel. Dia- 



gram (a) corresponds to an insertion of the momentum-dependent operator with coefficient C2 
diagrams (b)-(d) correspond to dressed OPGBE, while diagram (e) denotes an insertion of the 
m^-dependent operator with coefficient D2 ° . 



dressed by the LO amplitude. Also, there are two contributions from a single insertion of 
mq, with coefficients -D2A° ^23°^ which are also dressed by the LO amplitude. 

The scattering amplitude at NLO, ^i^^'"^'', is the sum of the contributions shown in 
Figs. and 
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It is straightforward to show that the individual contributions are 
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where d!^^"^ = D^^a^ + and /i is the renormalization scale. The PDS subtraction 



procedure has been used in defining the power-law divergent loop diagrams. 
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The partially-quenched amplitude in the ^5*0 channel can be found straightforwardly 
from the QCD amplitude by taking derivatives with respect to of the OPGBE contri- 
butions (see eq. ([T0|)) and by constructing the local operators that can contribute in the 
SU{4:\2)l ® SU{4:\2)ji EFT. The additional diagrams that contribute are shown in Fig. ^, 
and the scattering amplitude in the ^5*0 channel at NLO is 
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(a) 



(b) 





(d) 

FIG. 4. Additional contributions to the scattering amplitude in the ^5*0 channel at NLO, C(Q°), 
in the partially-quenched EFT. Diagram (a) denotes an insertion of the operator with coefficient 
-^2_B° ' diagrams (b)-(d) denote the dressed ry-exchanges. 
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where we have worked to LO in the relation between the quark masses and the meson masses. 
The explicit renormalization-scale dependence of the amplitude due to the log contri- 
bution is exactly compensated by the renormalization-scale dependence of the coefficient 
^23 (a*) to yield a //-independent amplitude. 

As there is no mixing between different partial waves in the ^5*0 channel it is straight- 
forward to determine the phase-shift, Sis^ = 5o,oo,05 from the scattering amphtude without 
approximation using eq. (p!3D, and, in turn, to construct the effective range expansion 
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Here a^^^°^ and r^^^°^ are the scattering length and effective range in the ^5*0 channel, respec- 
tively. The scattering length in the partially-quenched EFT is found to be 
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where 7 is a /x-independent linear combination of Cq ° and fj. that enters at LO in the 
expansion and must be determined from data. Furthermore, the effective range is found to 
be 
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In QCD with KSW power-counting, the scattering amplitude in the ^5*0 channel has 

and it has been found that the expansion is convergent. 



been determined up to NNLO |22 



However, the chiral expansion of the effective range parameters in this channel suggests 
that the convergence of the expansion is quite slow |31[]. Consequently, in order to have 
confidence in the chiral extrapolation of the partially-quenched amplitude and effective range 
parameters we have presented here, the NNLO amplitude (and even higher orders) should 
be computed in order to understand the convergence properties of the chiral expansion. 



C. The - '^Di Coupled Channels 

Due to the non-perturbative nature of OPGBE — particularly the chiral limit of the 
tensor force — in the ^Si — ^Di channel, the method for computation in this channel is 
fundamentally different from that in the ^5*0 channel and the higher partial waves where 
OPGBE can be included in perturbation theory. The details of the calculation of scattering 
lengths, phase shifts and bound state energies in the ^Si — ^Di coupled channels in QCD can 
be found in Refs. p3|,B2|, and we do not repeat them here. In the ^Si — ^Di coupled-channels, 



OPGBE generates both central and tensor potentials. 
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where a^, = (7^(1 — 2mldis/ gA) / {Sir f ) . The constant rfig is somewhat uncertain 
different extractions yielding -0.78 ± 0.27, -0.83 ± 0.06, -1.4 ± 0.24 Q and 
0.45 GeV^^ ||S^. Since the chiral limit of the potentials in eq. ( P^ contribute at LO 
shown in Fig. |], the m^-dependence of qa, f and are required at NLO 
of these observables has been studied extensively, the results of which can be found 
Refs. |3|J3|J33J3|, and up to NNLO it is found that 
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FIG. 5. Lippmann-Schwinger equation for the LO contribution to the scattering amphtude 
(large solid rectangle) in the ^Si — ^Di coupled-channels. The small solid circles denote an insertion 
of Cq or g^. The "o" appearing below the OPGBE diagram implies the chiral limit. 



/ = / 



(0) 



M 



N 



9A 



M 



(0) 
N 



(0) 
9 A 



1 - 



1 

~ 47r2(/(0))2 



ml log 



(PHYS) 



+ 



ml 



47r2(/(0)) 



m^ log 



^vr2(/(o))i 



i9T?ml 



+ 



Am 



9^' 



16 



(23) 



where mj^^^^^^ = 139 MeV, h = 4: A ± 0.2 0,0], ci 1 GeV"' H are mg-independent 

constants (we have explicitly separated the logarithmic contribution from I4). A complete 
analysis by Fettes |3^ of the vrA^ sector provides three determinations of diQ, diQ = —0.91 ± 
0.74, -1.01 ± 0.72 and -1.76 ± 0.85 GeV"^ At NLO there is a contribution from the 
chiral limit of TPGBE and from an insertion of a momentum dependent (p^) operator with 

coefficient C2^^\ as shown in Fig. |^. At this order there are two contributions arising from 
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together into D2 . The TPGBE potential in coordinate space has been computed in 



Ref. [17,40|, and in the chiral limit is given by 
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FIG. 6. Chiral limit of the crossed TPGBE diagram, deviations from the chiral limit of OPGBE, 
and the C2 (large solid circle) and D2 (large solid square) operators, all of which contribute at NLO 
in the ^Si — '^Di coupled-channels. The "o" appearing below a diagram implies the chiral limit. 



The singular nature of the tensor force has so far precluded regularization of this channel 
using dimensional regularization, or any other mass-independent regulator. We therefore 
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use a spatial square- well of radius R [^,^,0 , where the potential outside the square well, 
Vr^)(r;m.) = Vr^\is 
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The energy and mg-dependent potential, Vg^'^^\r] 171^,, k'^) = Vg^'^^\ inside the square well 
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where V^^'"^^ = V (Sg^) + ^n(^Si) + P^^(3si) and where we have again used the LO 
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relation between nig and the pion mass. V (Ss^), V (Sg^) and V (is-^) are constant potentials 



corresponding to the renormalized local operators with coefficients Cq , D2 and C2 
in the ^Si — ^Di coupled-channels, respectively. An identification can be made between 



the coefficients of the local operators, Cj;^^'' and Dl'^^', and the constant potentials of the 



square-wells that enter into eq. 



V^(3sj) and V^(3sj). For example 
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It is important to recall that there is implicit mg-dependence in this potential arising from 
the chiral expansions of qa, Mj\f and /, in addition to the explicit dependence from the 



contribution. Defining the two-component wavefunction \E' to be 
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where u{r) is the S-wave wavefunction and w{r) is the D-wave wavefunction, the regulated 
Schrodinger equation is 
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At LO in the partially-quenched theory one has contributions to the NN potential from 
the momentum and m^-independent four-nucleon operator with coefficient Cq and from 
the chiral limit of OPGBE. In fact, as discussed earlier, it is consistent to retain the full ruq 
dependence of OPGBE potential at NLO, which in PQQCD leads to 
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where it is implicit that both and ckq are defined in the partially-quenched theory, whose 
chiral expansions differ from those of QCD. Given that only the tensor component of OPGBE 
contributes at LO, and contributes only at NLO, only the partially-quenched expansion 
of is required for this NLO calculation. While one can straightforwardly construct all 
the operators that contribute to the process in PQxPT as one does in xPT, this is a te- 
dious procedure. Ultimately one arrives at the following relations appropriate for an NLO 
calculation 

g\ I 2mldis 2{mj.s - ml)disB Y gl 
= 1 , ao = ^-Fi , (32) 

where disB is an additional coefficient that must be determined from lattice calculations. 
In the chiral expansion of the OPGBE potentials in eq. (|3l| ) it is important to note that 
there is no contribution of the form The appearance of such a term would destroy the 
renormalization program we have constructed for QCD and it is encouraging that such a 
term does not arise from OPGBE in PQQCD either. 

The chiral expansions of /, M^r and qa are also required, and these are known. The 
chiral expansion of / at NLO is |^ 



/ = / 



(0) 



m|y , (msv\ , , 2 , , 2 



(33) 



where msv is the mass of a meson composed of one valence and one sea quark, at LO 
in the chiral expansion, and li^2 are coefficients that need to be determined from lattice 
calculations and are directly related to the constant I4 in QCD. At NLO, the nucleon mass 
receives contributions from counterterms only 

Mm = Mo + ciml + C2mls + ... , (34) 

where the ci^2 are to be determined from lattice calculations. The matrix element of the 
axial current is somewhat more complicated, as it depends upon how one extends the axial 
currents from QCD to PQQCD, as discussed in Refs. |4^J|,|l3]. For vanishing ghost and 
sea-quark axial charge y'^^\ the axial matrix element is found to be 

9A = + ^ (r/(°) - gfw^ ) + c(°) , (35) 



where 



(0)n3 (0) (0) 

/I X (0) (0) 

- = {9'^? (^^. + ^Lsv) - '-^-^ (5,r - gf) {K - Lsv) 

g{0) ^ ^2 _|_ ^2^ ^ ^gg-J 

and Lt, = ml- log (m^//i^), and Lsv = "^sv ^^S (^Iv//^^)- Generalization to arbitrary ghost 
and sea-quark axial charges is straightforward p3|,p!0|. The rj are constants that must be 



determined from the lattice. Unlike in the single nucleon sector |jTO[, we have not included 
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the A resonance as an explicit degree of freedom, as the A-nucleon mass sphtting is approx- 
imately the same as Atvtv, the scale at which the EFT is expected to break-down. If the 
range of the EFT is somehow extended to higher momenta (or if this one is shown to be 
valid at higher momenta that we presently expect), the A should be included in the theory 
explicitly, as the AA^vr intermediate states — in contrast with the AA states — are likely 
to make a sizable contribution that cannot be described by local counterterms at momenta 
beyond A at at- 



X 



FIG. 7. Additional contributions to the scattering amplitude in the ^Si channel at NLO, 0{Q^), 
in the partially-quenched EFT. The large solid square denotes an insertion of an m^-dependent 
operator with coefficient D^^^ ■ 



At NLO in the partially-quenched EFT there are potentially additional contributions 
from TPGBE. However, it is only the chiral limit of TPGBE that contributes at NLO and 
the chiral limit of PQQCD is the same as the chiral limit of QCD, by construction, and thus 
the TPGBE potentials are those given in eq. 0. Thus, the additional contributions at 
NLO have the same form as in the ^5*0 channel: a single insertion of iriq and the exchange 
of a single r], as shown in Fig. 0, which generates the potentials given in eq (PH) . The 
additional contribution from a single insertion of leads to a modification of the short- 



distance potential in eq. (pT]). The short distance potential in the partially-quenched theory 
is 



This completes the construction of a partially-quenched EFT describing the ^5*1 — ^Di 
coupled channels. The potential defined above is inserted into the Schrodinger equation 
in eq. (^) to generate observables. Due to the lack of partially-quenched lattice data in 
this channel we have not generated phase-shifts or scattering lengths. However, this is a 



straightforward procedure and requires only limited numerical work p3|j32 |. 



III. DISCUSSION 

In this work we have formulated the effective field theory required to describe the low- 
energy behavior of partially-quenched QCD in the two-nucleon sector. In fact, it is quite 
simple to construct the partially-quenched effective field theory from the known QCD results 



^In W counting one keeps the full TPGBE contribution at NLO |17], for which partial-quenching 
is somewhat more complicated 
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100 150 200 
m [MeV] 



300 



FIG. 8. The scattering length in the ^5i-channel in QCD as a function of the pion mass for char- 
acteristic strong-interaction parameters. The two shaded regions correspond to different allowed 
ranges oiD'^'^" that are both consistent with naive dimensional analysis. For a detailed discussion, 
see Ref. At the physical value of the pion mass the scattering length is a^^^^^ ~ +5.425 fm. 



and it is gratifying to see that one can obtain analytic results for many observables in the 
^5*0 channel and in the higher partial waves. While a numerical solution is required in the 
^Si — ^Di coupled channels, it amounts to a simple problem in non-relativistic quantum 
mechanics. 

One should be concerned about the range of sea and valence quark masses for which this 
theory converges. In QCD it is found that the NN EFT converges for and momenta less 
than of order A^ttv ~ 300 MeV, and one suspects that the same radius of convergence will 
exist in the partially-quenched theory. If this is indeed the case, lattice calculations will be 
required with meson masses of less than ~ 300 MeV in order to match to the EFT and use 
it to make predictions about nature. This is somewhat more restrictive than in the meson 
and single nucleon sectors and therefore one would like to see convergent results in those 
sectors before being confident in results obtained in the multi-nucleon sectors. 

It is very encouraging to see that partially-quenched calculations of quantities in PQxPT 
describing the dynamics of the PGB's are presently being performed [Q, and linear combi- 
nations of the Gasser-Leutwyler coefficients appearing at in xPT are being determined 
as a result. The situation is far less advanced in the two-nucleon sector. At present, NN scat- 
tering lengths must be extracted at finite volume using Liischer's formula, which expresses 
the energy of a two-particle state as a perturbative expansion in the scattering length divided 
by the size of the box [Q. There are (at least) two potential problems with this approach. 
First, one may worry that lack of unitarity in PQQCD may invalidate Liischer's formula for 
the NN scattering lengths. However, Liischer's formula is easily obtained in EFT(7f) |^ 



and one can convince oneself using the arguments of Ref. ||2^ that, while the pionful NN 



EFT described above is not unitary in PQQCD, EFT(7f) is unitary in PQQCD; all the effects 
of partial-quenching in EFT(7f) are in the coefficients of the contact operators. It follows by 
continuity that the NN scattering lengths can be extracted in a lattice simulation of PQQCD 
by using Liischer's formula and extrapolating to the physical quark masses using the for- 
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malism presented in this paper. A second worry is truly cause for concern: the S-wave NN 
scattering lengths are extremely large (which is understood as proximity to an infrared fixed 



) as compared to the sizes of state-of-the-art lattices: a*- ~ —23.714 fm and 
+5.425 fm. A recent study |^ of the pion- mass dependence of the NN scattering 
lengths in QCD suggest that the ^5*1 scattering length relaxes to natural values of ~ 1 fm as 



point 1^ 




■^i-.-'-r. ' 



/i - '-'t. -■ . .- . 1.-.'- -. ^ 



"0 50 100 150 200 250 300 

m [MeV] 



FIG. 9. The scattering length in the ^Sg-channel in QCD as a function of the pion mass 
for characteristic strong-interaction parameters. The two shaded regions correspond to differ- 



ed) 



ent allowed ranges of D2 
detailed discussion, see Ref. 



that are both consistent with naive dimensional analysis. For a 
2|. At the physical value of the pion mass the scattering length is 



-23.714 ± 0.013 fm. 



the pion mass is increased beyond ~ 200 MeV (see Fig. |^). One anticipates similar behavior 
in the partially-quenched theory. Given current uncertainties in strong interaction parame- 
ters, particularly -02^°'', it is at present unclear whether the same is true in the ^5*0 channel 
(see Fig. I). 

To our knowledge, a single lattice determination of the ^5*0 and ^Si NN scattering lengths 
in QQCD exists PH] at a pion mass of ~ 500 MeV. This pion mass is beyond the range 
of applicability of the EFT described in this paper, or of an analogous EFT that one can 
easily construct to describe QQCD. Given the unknown D2 coefficients that encode the short- 
distance quark-mass dependence in the S-wave channels one may question the motivation for 
an improved, partially-quenched simulation of the NN scattering lengths with pion masses 
within the NN EFT, since such simulation will, at best, simply determine the D2 operators. 
It has recently been shown ||23|j3^j49[] that to leading order in the NN EFT, and assuming 
perfect knowledge of the single-nucleon sector, the D2 operators determine the quark-mass 
dependence of the deuteron and, more generally, of dinucleon binding. Since small changes 
in ruq can in principle lead to drastic changes in the positions of nuclear energy levels, much 
attention has been given to light-element abundances predicted by big-bang nucleonsynthesis 
and to the abundance of isotopes produced by the Oklo "natural reactor" in the hope 
that these abundances can be used to constrain high-energy physics [|0|. This is powerful 
motivation indeed for improved lattice simulations of NN scattering lengths. We look forward 
to a significant lattice effort in the multi-nucleon sector. 
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